Anomalous exponents at the onset of an instability 
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Critical exponents are calculated exactly at the onset of an instability, using asymptotic expansion- 
techniques. When the unstable mode is subject to multiplicative noise whose spectrum at zero 
frequency vanishes, we show that the critical behavior can be anomalous, i.e. the mode amplitude 
X scales with departure from onset fi as (X) oc with an exponent P different from its deterministic 
value. This behavior is observed in a direct numerical simulation of the dynamo instability and our 
results provide a possible explanation to recent experimental observations. 



In the vicinity of a continuous phase transition, the 
amphtude of the order parameter, say M, increases with 
the departure from the critical point, say e, as a power 
law, i.e. M (X . Mean-field theories predict simple 
rational numbers for the exponent /3 (for instance 1/2 
for systems with cubic nonlinearities). It has been real- 
ized for a long time that, because of thermal fluctuations, 
the power law may differ from this mean-field prediction 
[H. The exponents are then said to be anomalous. Us- 
ing renormalization-group techniques, their value can be 
calculated as a perturbative expansion in the critical di- 
mension minus the spatial dimension of the system Q ■ 

Similarly, in the vicinity of a continuous instability in 
an out of equilibrium system, the amplitude of the un- 
stable mode, say X, grows with the departure from on- 
set, say as a power law {X) cx /i^ (where the angular 
brackets denote time-average). Dynamical systems ob- 
tained using normal form theory [3| provide simple ratio- 
nal values for /3 (usually 1/2 when the problem has the 
X — > —X symmetry, 1/4 at the tricritical point where 
the cubic nonlinearity vanishes and so on). Guided by 
the phase transition observations, one may expect that 
fluctuations shift the exponent /3 away from its mean-field 
value. Somehow surprisingly, the overwhelming major- 
ity of experiments on instabilities reports simple rational 
values in agreement with the mean-field prediction for 
anomalous exponents seem not to be measured in this 
context 0, 0|- In a recent experiment in a turbulent 
flow of liquid sodium, the dynamo instability has been 
observed and some measurements indicate that the first 
moment of the magnetic field displays an exponent 0.77 
. It is possible that experimental biases arc responsible 
for this observation: the instability is slightly imperfect 
and the numerical value of the exponent is then highly 
sensitive to the accuracy of determination of the onset. 
Another appealing possible explanation is that the tur- 
bulent fluctuations of the flow lead to the anomalous ex- 
ponent With the latter in mind, we now describe a 
canonical model that leads to anomalous behavior similar 
to the one measured in the dynamo instability. 

In the dynamo context, the turbulent fluctuations act 
as a multiplicative term in the equation for the magnetic 
field. In contrast to the case of equilibrium phase transi- 



tion where additive thermal fluctuations prohibit phase 
transition in small dimensions, bifurcations are not de- 
stroyed by multiplicative fluctuations even for small (pos- 
sibly zero) dimensions. We thus start with a zero dimen- 
sional system subject to multiplicative noise. For a mul- 
tiplicative white noise, on-off intermittency is a generic 
behavior close to the threshold of instability [1, 0| . Then 
the averaged amplitude scales as (X) ex. ^. Although 
the exponent differs from the mean-field prediction, its 
value /3 = 1 is in disagreement with the one measured in 
the experimental dynamo. It has been shown that on-off 
behavior is observed when the departure from onset is 
smaller than half of the value of the noise spectrum at 
zero frequency [lol |. In the dynamo experiment, on-off 
intermittency is not observed. We suggest that it is due 
to the absence of noise component at zero frequency and 
to strengthen this hypothesis, we investigate the effect of 
a noise whose spectrum at zero frequency vanishes. We 
thus consider the dynamics of the unstable mode X given 

by 



X = fix - a:"+i 



YX, 



(1) 



Here C is a Gaussian white noise with {C{t)({t')) = 
25{t — t'). is a (potential) function of Y and the 
subindex denotes differentiation with respect to this vari- 
able. Y acts as a multiplicative noise (for X) whose 
frequency-spectrum is controlled by the function F{Y). 
When the potential F is such that the second moment of 
Y is finite, the spectrum of Y vanishes at low frequency 
(it behaves as the square of the frequency /, for small 
/). Standard estimates of the effect of noise on the onset 
of instability (for instance by calculating the evolution of 
the ensemble average of log X from the linear part of the 
first equation [lH) show that the onset of instability of 
the solution X = is not affected by the noise and re- 
mains at /i = 0. In contrast, the non-linear regime above 
onset is strongly affected. We display in fig. [1] time se- 
ries of X for different functions F in the vicinity of the 
onset of instability (unless otherwise stated, numerical 
simulations are performed in the case of cubic nonlinear- 
ities: n = 2). For panel (a) we used white noise, F = 
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FIG. 1. Time series of the solution of eq. ([ij for ^ — 0.01. Top: linear scale and bottom in log scale, a) F = corresponding 
to a white noise; 

b) F^u = 7^^/2, Ornstein-Uhlenbeck noise with 7 = 1.5; c) -F^jy = i^\Y\, with v = 0.75. Note the differences 
in the y-coordinate values. 



and on-off intermittency is observed: short bursts of fi- 
nite amplitude (on-phases) alternate with long durations 
with negligible amplitude (off-phases). In panel (b) the 
case F = F^^ = is presented. For this choice Y 

is the Ornstein-Uhlenbeck process. There is no off-phase 
and we expect a behavior for the moments that differs 
from the one of on-off intermittency. Panel (c) displays 
a time series for F = F^^^ = I'lY], that results in an 
intermediate behavior. 

In fig. [2] the first moment is displayed as a function 
of ^ for the two functions F^^ with 7 = 0.2 and 
with various values of v. For the F^^ case we observe for 
/i e [3.10^^,10^^] an evolution that seems compatible 
with a power law. A best fit determination of the asso- 
ciated exponent results in the value 0.69, thus different 
from 1 and 1/2. However when ^ is very small, the slope 
changes and the deterministic exponent 1/2 is recovered: 
the apparent anomalous behavior disappears at critical- 
ity |12| . This is confirmed by a perturbative expan- 
sion performed on the Fokker-Planck equation (not pre- 
sented here). This expansion predicts that X is concen- 
trated around the value X* at which a weighted average 
of the non-linear effect balances the linear growth rate 
^ = X* " U{Y) exp (nY) dY, where n(r) oc e'^ is 
the stationary probability density of Y. Thus, in this 
case and for n = 2, the first moment scales as as 
observed numerically. 

A simple potential F for which this expansion can 
break down is F^^ = i^lYl. Indeed, if > n the ex- 
pansion holds resulting in normal scaling /3 ~ l/n but 
breaks down (because X* vanishes) when i' < n. In fig. 
[2I where the first moment for this potential is displayed, 
we observe that {X) oc jjL for small v and {X) oc y/Ji 
for large v. Anomalous behavior with exponent between 
1/2 and 1 is observed for v of order 1. In this regime 
and in contrast to the case, the exponent remains 




FIG. 2. First moment (X) as a function of /i for the solution 
of eq. O for F = Fan with (■) v = 0.125, (♦) v = 1.125, (•) 
V — 1.375, (a) u — 1.75, (ir) v = 2.5. The data are presented 
in loglog scale and have been normalized by their value at 
= 0.01. The results for F = Fou with 7 = 0.2 (o) are 
presented and shifted for comparison. The thick continuous 
lines indicate the exponents 1/2, 0.69 and 1. 



anomalous for the smallest achievable values of /i. This 
numerical result is confirmed by a new perturbative ex- 
pansion that we now sum up. 

Using Vt ~ log AT — y — log/i/n, the Fokker-Planck 
equation for P the stationary probability density function 
(p.d.f.) of VL and Y is 

= -^l^n[l - e"^^+"^)P + dYiFyP) + d^P . (2) 

Since the derivative in VL is multiplied by a small pa- 
rameter (we are interested in the limit /i — > 0), we in- 
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troduce a WKB-like expansion and search for P(r2, y) = 
exp [X]m=-i A*™'^"*] J where the first term S-i depends 
only on f2. At lowest order we obtain 



a2ro + 9y(Fyro) + 5-i,o(e 



nfl+nY 



l)ro=0, (3) 



where vq = exp[5'o]. This equation can be solved exactly 
for positive and negative Y. The two solutions are then 
matched at y = which selects the value of S-i 



(4) 



where A 



-45_i,o/n2 



jv? — and and 



if„ are modified Bessel functions of order k. The solution 
for tq is then 



(5) 

The amplitude ^(51) is determined from the solvability 
condition at next order. Up to this order, we have then 
obtained the expression P = exp[/z~^S'_i(r2)]ro(r2, F) 
where all the dependence in is in the exponential. As 
displayed in fig. [3l this asymptotic result is in good agree- 
ment with the numerical simulations of the Langevin 
equations ([T]). 
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FIG. 3. Probability density function 11 of + log(/i)/n = 
log(X) — Y . The continuous line is the theoretical prediction 
and the symbols are the PDF calculated from the numerical 
solutions of the Langevin equation. Here v = 1 and the three 
curves are associated (from left to right) to ^ = 1.7810^''', 
H = 3.16 10"'' and n = 5.6 10-^ 



From this formulation, we can calculate the moments. 
The exponential term acts as a cut-off for large VI and 
is of the form exp (— exp (n vVl/ (n — i^))). Therefore 
for jjL ^ Q and v < n , only very negative f2 have to be 
considered. In this hmit the amplitude A{^)Kf^ [Ae"^/^] 
tends to a constant and. after several standard estimates 



of the asymptotic behavior of the Bessel functions, we 
obtain for u < n 



/3 



(6) 



We tested our prediction by numerically calculating 
the first moment for different values of v and for n = 2 
and n — i. The results are shown in figure 01 For all 
cases the predictions are within the error-bars of the nu- 
merically calculated values of /?, and thus the predic- 
tions are verified. To discuss one particular value, the 
numerically computed exponent for v ~ and n = 2 is 
j3 = 0.66 ± 0.02 which is in perfect agreement with the 
theoretical prediction 2/3. We have also performed sev- 
eral numerical simulations using potentials of the form 
F = —vyjY^+Y^ . We have observed that only the 
behavior of F for large values of \Y\ is important. In 
other words, the universality classes of the problem (i.e. 
the models having the same critical exponents) are deter- 
mined by the behavior of the tails of n(y). Incidentally, 
this shows that the anomalous scaling is not caused by 
the non-analyticity of F at F = 0. 

At this stage, we emphasize that our perturbativc ex- 
pansion (in ^) allows to calculate an exact (non pertur- 
bativc) expression for the value of the anomalous expo- 
nent. This exponent transitions from its on-off value 1 
for V <\ to its deterministic value \/n for v > n. In the 
simple case of cubic nonlincarities, we predict an expo- 
nent between 1/2 and 1. Interestingly enough, the scaling 
reported in the dynamo experiment belongs to this range. 
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FIG. 4. Exponents of the first moment as a function of u 
for (□): n — 2 and (o): n = 3. The continuous lines are the 
theoretical predictions and the symbols are obtained from the 
numerical solutions of the Langevin equation. 



We have focused here on the first moment of the un- 
stable mode. The behavior of higher moments is also of 
interest. It can be characterized by the set of exponents 
j3p defined by {X^) oc . In the absence of fluctuations 
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or at usual equilibrium phase transitions, monoscaling 
is observed which means that /3p = pf3i. The situa- 
tion is richer here: there is no linear relation between 
the exponents (for instance it can be easily proved that 
/3„ = 1). Thus the solutions of model ^ display multi- 
scaling. This is related to the complex structure of the 
p.d.f. of X. In particular, it cannot be expressed as 
a simple one-parameter distribution characterized by its 
first moment in contrast to the scaling hypothesis close 
to the critical point of an equilibrium phase transition 

Another important issue is the effect of spatial dimen- 
sion. The model ([T]) is zero dimensional {X only depends 
on time and not on space) while the magnetic field in 
magnetohydrodynamics (MHD) depends on three spatial 
dimensions. Analytical predictions for the critical behav- 
ior at larger (non-zero) dimensions would be of great in- 
terest but are still out of reach at present. To investigate 
further the pertinence of our model to the dynamo insta- 
bility, we have performed direct numerical simulations 
of the MHD equations. To increase our control on the 
velocity temporal behavior, we used the infinite Prandtl 
number limit 15[. In this limit the velocity is slaved to 



an external mechanical forcing and the Lorentz force 



b - Vb- VP, 



where b is the magnetic field and F is the body force. 
It is proportional to the ABC fiow F = A„[5sin(z) + 
2 cos(j/) j_2 sin(a;) -I- 5 cos(2:), 2 sm{y) + 2 cos(a;)] (see for in- 
stance [l6[)- is an amplitude that changes every 
time interval r based on a discrete version of our model 
An+i ^Ao + (y„+i - y„) and r„+i = r„ - TF(r„) + r„ 
where r„ is a random number. The magnetic field satis- 
fies the induction equation 

8h 

— = V X (u X b) + iim-^V^b. 

ot 

The MHD equations were solved in a periodic box of 
size 27rL using a standard pseudo-spectral code [ij] on 
a grid 32'^. The magnetic Reynolds number defined by 
Rm = (II^IP)^^^^/^ was varied above the onset value 
Rrric — 11.65. In fig. [5] (a), we display time series of the 
magnetic energy and note that they are similar to those 
presented in fig. [1] The first moments are displayed in 
fig. [5] (b) for several values of ly. We observe that the 
exponent of the first moment decreases from 1 to 1/2 
when v increases. Estimates of the exponent are compu- 
tationally demanding so that a quantitative comparison 
with our model is out of reach. Nevertheless, the results 
reported here support the robustness of the behavior wc 
have identified. 

To summarize, we have presented a simple model that 
results in anomalous exponents which lie between the 
deterministic value and the on-off intermittent one. The 
exact value of these exponents was calculated using an 
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FIG. 5. (Top) Time series of the space-averaged magnetic 
energy Et = B"^ above the dynamo onset for (from top to 
bottom) F — Fou with ■y — 1, F = Fan with u = 0.4 and 
a white noise (see discussion in the text). The curves have 
been shifted for clarity. (Bottom) First moment as a function 
of the departure from onset {Rm — Rrric) / Rrric for F = Fan 
with (A): u ^0, (+): u = 0.1, (□): u = 0.4 and (o): u = 0.8. 
The two thick lines indicate the exponents 1/2 and 1. 



asymptotic expansion. The model emphasizes the role 
of the noise spectrum at zero frequency. It remains to 
be understood whether and when turbulent fluctuations 
can be modeled as the noise considered here 



ITf. In addi- 



tion, how such a noise affects other phase transitions and 
whether the present expansion can capture other critical 
exponents are interesting open questions. 
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